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ABSTRACT 

The  instantaneous  force  acting  on  a  circular  cylinder 
immersed  in  a  unidirectional  periodic  flow  characterized  by 
U  =  V  +  V.  sin  2iTft  was  investigated.   The  total  force  was 
divided  into  three  components  consisting  of  a  steady-state 
drag;  an  oscillating  drag  which  is  a  function  of  the 
frequency  parameter   X  =  fd/V  ;   and  an  inertial  force  v;hich 
is  also  a  function  of  the  frequency  parameter.   The  three 
resistance  coefficients  associated  with  the  force  components, 
Cj^(O),  C^(X),  and  Cj^(X)  were  calculated. 

It  was  found  that  the  resistance  coefficients  are  not 
materially  affected  by  the  free-stream  oscillations  whenever 
the  frequency  of  the  oscillation  is  lower  than  the  natural 
frequency  of  the  vortex  shedding.   However,  when  X  approaches 
the  Strouhal  frequency  of  vortex  shedding,  the  inertia 
coefficient  Cj^  is  significantly  reduced,  and  likewise  when 
X  approaches  the  value  of  twice  the  Strouhal  number,  C^(X) 
is  markedly  reduced.   In  regions  below  this  range,  the 
inertial  force  becomes  negligible  and  the  total  resistance 
could  be  represented  accurately  by  the  coefficients  ^^,(0) 
and  C_^(X)  with  Cj^  fixed  at  its  theoretical  unseparated-flow 
value,  2.0. 
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NOMENCLATURE 

Symbol 

A  effective  cylinder  flow-blockage  area  (A  =  6.1) 

a  instantaneous  acceleration 

A  cylinder  cross-sectional  area  (A  =  Trd  /4) 

A^„  cross-sectional  area  of  test  section  (Arpc  =  h£) 

Cp.  steady  state  drag  coefficient 

C_^(0)  coefficient  for  steady  state  portion  of  drag 

C_(X)  coefficient  for  oscillating  portion  of  drag 

C^(X)  coefficient  for  inertia  force 

d  diameter  of  test  cylinder 

e  tunnel  correction  factor 

F  instantaneous  total  force  acting  on  the  test 
cylinder 

f  frequency  of  oscillation 

P^  drag  forces  acting  on  the  test  cylinder 

F-j.  inertia  forces  acting  on  the  test  cylinder 

h  height  of- tunnel  test  section 

K^  constant  (K,  =  '-^pV^A) 

K2  constant  (K^  =  VV^Ap) 

K-  constant  (K   =  hV^^kp) 

K^  constant  (K^  =  2TrfpA  iV   ) 

I  width  of  tunnel  test  section;  test  cylinder  length 

n  frequency  of  vortex  shedding 

P  pressure 

p  height  of  dynamic  pressure  curve  at  t  =  T/4 


q  measured  velocity 

R  test  cylinder  radius 

r  height  of  dynamic  pressure  curve  at  t  =  3T/4 

Re  Reynolds  number  (Re  =  Ud/v) 

Re  Reynolds  number  (Re  =  Vd/v) 

St  Strouhal  number  (St  =  nd/V) 

T  period  of  oscillation 

t  time 

U  instantaneous   velocity    (U  =   V  +  V.    sin2iTft) 

V  free  stream  velocity 

V.  amplitude  of  velocity  oscillation 

V  average  velocity 

X  horizontal  distance  between  pressure  taps  1  and  3 

V  horizontal  distance  between  pressure  tap  3  and 
center  of  test  cylinder 

X  normalized  frequency  of  oscillation  (X  =  fd/V) 

V  fluid  kinematic  viscosity 
p  fluid  density 


I.   INTRODUCTION 

A.   EXPERIMENTAL  JUSTIFICATION 

The  subject  of  forces  acting  on  blufr  bodies  immersed 
in  time-dependent  flows  has  been  and  will  continue  to  be  of 
interest  to  fluid  dynamicists,  aerodynamicists,  and  practicing 
engineers  for  the  special  reason  that  in  nature  neither  the 
body  nor  the  fluid  which  surrounds  the  body  is  ever  in  a 
state  of  steady  motion.   Even  the  flow  behind  a  bluff  body 
moving  steadily  through  a  fluid  is  accompanied  by  large 
scale  unsteadiness.   Thus,  any  type  of  unsteadiness  of  the 
ambient  flov;  and/or  the  motion  of  the  body  introduces 
additional  changes  in  the  characteristics  of  the  flow  and 
its  analysis. 

The  drag  and  inertial  forces  are  interdependent  as  well 
as  time  dependent,  and  the  resistance  coefficients  obtained 
in  unseparated  flows  are  not  applicable  to  separated  cases. 
Although  indirect,  the  role  of  viscosity  is  paramount  in 
that  its  consequences  are  separation,  vortex  formation  and 
shedding,  and  resultant  alternations  in  the  virtual  mass. 
It  is  thus  clear  that  it  is  necessary  to  determine  the 
relationships  between  various  resistance  components  in  terms 
of  the  unsteadiness  of  the  ambient  flow  and/or  body  motion, 
geometry  of  the  body,  the  degree  of  upstream  turbulence, 
past  history  of  the  flow,  etc.  [Ref.  1]. 


However,  the  understanding  of  the  behavior  of  time- 
dependent  flows  about  bluff  bodies  may  be  advanced  only  by 
considering  the  relatively  more  manageable  cases  and  gradually 
integrating  the  information  so  obtained.   This  led  to  the 
present  experimental  investigation  of  the  unidirectional 
periodic  flow  about  circular  cylinders. 

B.   SURVEY  OF  THE  PREVIOUS  INVESTIGATIONS'^ 

A  review  of  some  of  the  previous  investigations  on  time- 
dependent  flows  may  best  be  presented  by  dividing  the  motions 
under  consideration  into  various  but  admittedly  arbitrary 
categories. 

1.   Theoretical  Analysis  of  Unseparated  Time-Dependent 
Laminar  Flows 

The  study  of  the  theory  of  separation-free  time- 
dependent  laminar  flows  has  enjoyed  particular  attention, 
partly  due  to  its  practical  significance  and  partly  due  to 
its  relative  mathematical  simplicity.   Surveys  by  Stewartson 
[Ref.  2],  Stuart  [Ref.  3],  and  Rott  [Ref.  4]  reflect  the 
current  level  of  understanding. 

The  unseparated  class  of  unsteady  f lov7s ,  most  of 
which  results  from  the  unidirectional  or  periodic  acceleration 
of  bodies  in  an  infinite  or  bounded  fluid  medium,  gives  rise 
to  an  induced  mass  which  must  be  added  to  the  real  mass  of 
the  body.   In  general,  the  virtual  inertia  tensor  depends 
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on  the  shape  of  the  body,  the  nature  of  the  fluid  medium, 

the  orientation  of  the  moving  object  with  respect  to  the 

direction  of  motion,  and  the  depth  of  relative  submergence 

from  a  free  surface  and/or  solid  boundary. 

The  effects  of  viscosity  and  separation  on  added 

mass  or  the  dependence  of  added  mass  of  a  given  body  on  the 

characteristics  of  the  time-dependent  separated  flow  cannot 

yet  be  theoretically  evaluated  nor  can  it  be  experim.entally 

separated  from  the  total  resistance  experienced  by  the  body. 

2.   Small  Amplitude  "Vibration  of  Bodies  in  a  Liquid 
Othervjise  at  Rest  and  the  Initial  Instants  of 
Acceleration  of  Bodies  from  a  State  of  Uniform 
Velocity 

Vibratory  motion  has  often  been  used  [Refs.  5-8]  as 
one  of  the  experimental  means  to  determine  the  virtual  mass 
coefficients.   The  technique  has  been  restricted  to  amplitudes 
of  motion  which  are  so  small  that  either  separation  does  not 
occur  or  the  characteristics  of  separation  do  not  change 
due  to  the  imposed,  high-frequency  vibration.   Thus,  the 
results  obtained  are  not  applicable  to  occurrences  in  which 
the  duration  of  unsteady  flow  in  one  direction  is  long  enough 
for  separation  to  occur. 

It  is  not  necessary  that  the  body  be  initially  at 
rest  to  determine  the  added  mass  provided  that  one  starts 
with  a  steady  ambient  flow  about  the  body  or  with  the  body 
moving  at  constant  velocity  in  a  fluid  othervjise  at  rest. 
Hamilton  and  Lindell  [Ref.  9]  have  shown  that  the  added-mass 
coefficient  of  a  sphere  determined  by  imposing  a  uniform 
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acceleration  onto  its  steady  state  motion  is  equal  to  its 
theoretically  determined  value  as  well  as  to  its  experimen- 
tally determined  value  obtained  by  small-amplitude,  high- 
frequency  oscillations. 

Rayleigh  [Ref.  10]  has  shown  that  if  the  acceleration 
continues,  i.e.  V  =  V  +  kt  ,  both  the  drag  and  the  added-mass 
coefficients  change  with  time. 

The  foregoing  arguments  suggest  that  the  force 
exerted  by  the  fluid  on  a  body  immersed  in  a  time-dependent 
flow  might  be  expressed  as 

F  =  hC^pk\v\Y   +  Cj^  dV/dt  +  history-dependent    (1) 
drag  and  inertia  forces 


in  which  A  is  the  projected  area  of  the  body  and  m„  the  mass 

of  the  fluid  displaced  by  the  body.   The  drag  coefficient  C^ 

is  to  be  taken  equal  to  the  quasi-steady  state  resistance 

coefficient  and  C„  to  be  taken  equal  to  that  measured  or 

calculated  from  the  potential  theory  when  the  value  of  the 

history-dependent  force  is  zero  [Ref.  9]. 

3.   Bodies  Subjected  to  Large  Amplitude  Harmonic 
Oscillations  in  a  Fluid  Otherv/lse  at  Rest 

A  number  of  significant  studies  [Refs.  11-15]  were 
made  to  determine  the  forces  acting  on  bluff  bodies  under- 
going harmonic,  large-amplitude  oscillations  in  a  liquid 
otherwise  at  rest. 
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Odar  and  Hamilton  [Ref.  11]  measured  the  force  acting 
on  a  sphere  subjected  to  harmonic  oscillations  and  proposed 
a  force  equation  comprised  of  three  parts:   a  steady  drag 
component  where  the  drag  coefficient  is  the  well-established 
coefficient  for  steady  translation  of  a  sphere;  an  inertial 

component  whose  coefficient  is  determined  experimentally  as 

2 
a  function  of  V  /ad  (a:  acceleration;  d:  diameter  of  sphere); 

and  a  history-dependent  component  whose  coefficient  also  is 

2 
a  function  of  V  /ad. 

^.   Bodies  Subjected  to  Wave  Motion  with  Zero  Mean  Velocity 

Considerable  work  [Refs.  16-21]  has  been  done  on  this 

type  of  time-dependent  flow  in  an  attempt  to  predict  the  wave 

forces  on  piles  and  other  submerged  structures.   In  general, 

the  total  force  acting  on  the  body  is  assumed  to  be  composed 

of  a  velocity-dependent  and  an  acceleration-dependent  force. 

The  resulting  equation,  known  as  the  Morison  equation  [Ref. 

16],  is  given  by 

F  =  i^Cj^pA|V|V  +  Cj^^  dV/dt  (2) 

where  C.  =  1  +  C,  C  being  the  added  mass  coefficient.   Experi- 
ments show  that  C^   and  C   depend  on  time  and  show  considerable 
scatter  about  their  assumed  mean  values  [Refs.  17,  20]  over 
a  cycle. 


12 


5.   Bodies  Subjected  to  Unidirectional  Acceleration  in 
a  Fluid  Othervjise  at  Rest  or  Unidirectional  Unsteady 
Flow  About  Bodies  Held  at  Rest 

Substantial  effort  has  been  made  to  determine  the 
components  of  force  acting  on  a  body  accelerating  unidirec- 
tlonally  in  a  fluid  otherwise  at  rest.   Iverson  and  Balent 
[Ref.  22]  and  Keim  [Ref.  23]  towed  spheres,  disks,  and  cylin- 
ders in  unsteady  motion  through  still  water.   It  was  concluded 

that  the  total  instantaneous  force  on  these  objects  could  be 

2 
described  by  a  single  coefficient  of  the  form  C  =  2F/pAV 

and  that  this  coefficient  was  a  function  only  of  a  so-called 

2 
acceleration  modulus,  ad/V  .   Keim,  hov/ever,  detected  an 

apparent  Reynolds  number  effect  for  cylinders. 

Laird  et.  al.  [Ref.  2^]  measured  the  forces  acting 
on  accelerating  cylinders  and  found  strong  evidence  of  devia- 
tion of  the  drag  coefficient  from  the  accepted  values  for 
uniform  motion.   They  further  found  that  the  acceleration 
modulus  did  not  correlate  the  resistance  coefficient  near 
boundary-layer  transition. 

Sarpkaya  and  Garrison  [Ref.  22]  working  with  uni- 
directional flow  with  constant  acceleration  about  cylinders 
and  plates  have  found  that  both  C-^  and  C^  depend  on  the 
history  of  motion,  that  the  Morison  equation  (Eqn.  2)  is 
valid  only  for  flows  with  constant  acceleration,  and  that 
the  acceleration  modulus  can  correlate  the  data  only  if  the 
acceleration  is  kept  constant. 
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6.  Forced  or  Self-Induced  Transverse  Oscillations  of 
a  Body  in  a  Fluid  in  Steady  Motion 

Experiments  have  been  directed  to  the  effects  that 
attend  the  oscillations  of  a  body  in  the  plane  of  the  lift 
force.   This  type  of  interaction  becomes  particularly  impor- 
tant when  the  oscillation  of  the  structure  is  the  necessary 
condition  for  the  generation  of  the  exciting  forces.   The 
most  pertinent  information  from  these  experiments  is  that 
the  lift  and  drag  forces  act  on  the  oscillating  cylinder  at 
the  vortex  shedding  frequency  and  twice  the  shedding  fre- 
quency, respectively,  provided  that  the  driving  frequency  on 
the  cylinder  is  appreciably  different  from  the  shedding  fre- 
quency.  When  the  forcing  frequency  of  the  cylinder  approaches 
the  shedding  frequency,  the  natural  shedding  frequency  is 
lost  and  it  "locks-in"  to  the  forcing  frequency.   This 
synchronization  persists  over  a  range  of  frequencies  which 
may  be  termed  the  "range  of  synchronization"  [Ref.  26]. 

7.  Unidirectional  Oscillatory  Flow  About  Bodies  at  Rest 
and  Steady  Uniform  Flov;  About  Bodies  Subjected  to 
Streamwise  Oscillations 

Relatively  few  studies  have  been  conducted  in  this 
category.   Chen  and  Ballengee  [Ref.  2?]  examined  the  vortex 
shedding  from  circular  cylinders  in  an  oscillating  free 
stream.   Their  results  suggested  that  "in  an  oscillatory 
freestream  of  3  HZ  and  Reynolds  number  up  to  40,000,  the 
vortex  shedding  from  a  circular  cylinder  responds  instantan- 
eously to  the  freestream  variations"  and  that  "the  instan- 
taneous Strouhal  number  stays  sensibly  constant  at  0.2  ±  0.01." 
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Hatfield  and  Morkovin  [Ref.  28]  studied  the  effect 
of  an  oscillating  free  stream  on  the  unsteady  pressure  on  a 
circular  cylinder.   They  have  found  that  there  is  no  signi- 
ficant coupling  between  the  small-amplitude  free-stream 
oscillations  and  the  vortex  shedding.   Their  results  would 
suggest  that  the  drag  coefficient  associated  with  the  mean 
flow  would  essentially  remain  constant  at  its  steady  state 
value.   On  the  other  hand,  Mercier  [Ref.  29]  who  subjected 
cylinders  to  large  streamwise  oscillations  found  that  the 
average  drag  coefficient  significantly  increases  with  fd/V 
and  that  the  rate  of  increase  depends  on  the  amplitude  to 
diameter  ratio.   Combined,  these  results  would  suggest  that 
the  degree  of  coupling  between  the  frequency  of  oscillation 
and  the  vortex-shedding  frequency  and  hence  the  forces  acting 
on  the  body  strongly  depend  on  the  amplitude  of  oscillation. 

Davenport  [Ref.  30]  subjected  bluff  bodies  to  small 
amplitude  oscillations  in  a  water  flume  and  evaluated  the 
drag  and  inertia  coefficients  through  the  measurement  of  the 
rate  of  damping  of  the  amplitude  of  oscillations  of  the 
bodies,  i.e.  v/ithout  measuring  the  forces  acting  on  the 
bodies.   Using  the  frequency  parameter  fd/V  to  correlate 
the  data,  he  found  large  variations  in  C„  and  C,,. 

C.   SCOPE  OF  THE  PRESENT  STUDY 

The  present  investigation  is  limited  to  an  experimental 
determination  of  the  instantaneous  forces  acting  on  a 
circular  cylinder  immersed  in  a  unidirectional  periodic 
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flow   characterized  by 


U  =   V  +   V^  sin   27Tft.  (3) 


Its  theoretical  aspects  are  structured  along  the  lines  of 
the  unpublished  work  by  Sarpkaya  [Ref.  31]  conducted  at 
Aerodynamische  Versuchsanstalt-Gottingen  with  respect  to 
the  supposition  that  the  instantaneous  force  acting  on  an 
object  immersed  in  a  periodic  flow  strongly  depends  on  the 
parameter  fd/V  and  to  varying  degrees  on  Vd/v  and  V,/V 
where  f  is  the  frequency  of  oscillation  and  d  the  cylinder 
diameter. 

The  range  of  these  parameters  during  the  experimental 
work  was : 

0.0001  <  fd/V  <  0.05 
11,000  <  Vd/v  <  63,000 
0.085   <  V^/V  <  0.69 

As  presented  in  Section  IV,  the  parameter  fd/V  is  extended 
to  a  value  of  1.5  using  the  data  of  others. 
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II.   APPARATUS  AND  PROCEDURE 

A.   APPARATUS 

1.  Water  Tunnel  and  Oscillating  Mechanism 

The  experiments  were  conducted  in  a  500-gallon- 
capacity,  recirculating-type  water  tunnel  as  shown  in  Figure 
1.   V/ater  was  circulated  by  a  low-rpm,  high-capacity,  centrif- 
ugal pump  which  was  capable  of  providing  test-section  mean 
velocities  of  2-25  feet  per  second.   The  test  section 
consisted  of  a  rectangular  conduit  4  inches  wide,  8  inches 
high,  and  l6  inches  long. 

Periodic  variations  in  water  velocity  were  provided 
by  a  rotating  butterfly  valve  driven  by  a  variable  speed 
motor  capable  of  a  range  of  0-^0  rpm  (Fig.  2).   The  magni- 
tude of  the  velocity  fluctuations  was  dependent  upon  the 
loss  in  effective  head  through  the  valve.   As  this  is  partially 
a  function  of  the  valve-to-conduit  area  ratio  [Refs.  35,36], 
an  access  plate  was  placed  on  the  tunnel  in  order  to  change 
the  effective  valve  area  by  adding  concentric  rings  to  the 
valve  perimeter  without  removing  or  changing  the  valve  body. 

Associated  with  the  tunnel  is  a  15-gallon  head  tank 
located  12  feet  above  the  test  section  centerline. 

2.  Test  Specimens 

The  test  specimens  consisted  of  polished  aluminum 
cylinders  placed  across  the  width  of  the  test  section  perpen- 
dicular to  the  direction  of  flow  (Fig.  3) •   A  small  bore 
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hole  drilled  radially  from  the  cylinder  surface  at  the  mid- 
section served  as  the  total  pressure  tap.  The  test  cylinder 
was  mounted  with  one  end  held  by  a  self-aligning  bearing  fit 
in  one  plexiglass  window,  with  the  other  end  passing  through 
the  other  window  to  a  cantilever  beam  force-measurement 
system.  Two  cylinders  of  diameters  0.75  and  1.0  inches  were 
used. 

3.   Instrumentation 

The  total  force  acting  on  the  cylinder  within  the 
test  section  was  measured  directly  by  a  strain-gage  fitted 
cantilever  beam  (Fig.  k) .      The  bridge  circuit  consisted  of 
four  active  strain  gages,  with  the  output  signal  being  fed 
to  a  dual-channel  strip  chart  recorder. 

Three  pressure  taps  were  utilized  for  velocity  and 
acceleration  calculations.   Two  variable-sensitivity, 
differential-pressure  transducers  connected  to  two  dual- 
channel  strip  chart  recorders  were  used  for  this  purpose. 
Each  transducer/recorder  system  was  calibrated  with  a  water- 
manometer. 

B.   PROCEDURE 

1.   Measurem.ent  of  Force 

A  continuous  trace  of  the  total,  instantaneous  force 
acting  on  the  cylinder  was  obtained  from  the  strain-gage 
fitted  cantilever  beam  (Fig.  4)  which  was  wired  to  a  strip 
chart  recorder.   This  system  was  calibrated  by  static  loading 
of  the  cylinder.   Three  different  cantilever  beams  were  used 
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during  the  course  of  the  investigation,  each  beam  being 
several  times  more  sensitive  than  the  previous  one,  in  order 
to  increase  the  range  of  experimentation. 

2.  Measurement  of  Frequencies 

Each  strip  chart  recorder  has  the  capability  of 
variable,  controlled  paper  speeds.   Thus,  with  chart  speed 
and  distance  between  corresponding  points  on  the  trace 
known,  the  period  of  the  cyclic  signal  could  be  determined. 
The  frequency  of  oscillation  was  then  simply  the  reciprocal 
of  the  period. 

3.  Measurement  of  Acceleration 

Acceleration,  with  the  exception  of  a  scale  factor, 
was  measured  directly  from  the  differential  pressure  taken 
from  two  static  pressure  taps  placed  a  knov/n  distance  apart. 
Neglecting  friction  losses,  the  equation  of  motion  yields 
(Fig.  5), 


(P^  -  P.3)A^2  =  ma  (4) 


where  m  is  the  mass  of  fluid.   Writing  m  =  p  X  Arp^  in  (^), 
one  has 


(P^  -  P3)A^g  =  p  X  A^g  a  (5) 


or 


a  =  AP^_3/pX.  (6) 
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^.   Measurement  of  Velocity 

2 
The  dynamic  pressure,  In  the  form  pU  /2,  was 

determined  using  the  differential  pressure  from  the  total 

and  static  pressure  taps  2  and  3,  respectively,  (Fig.  5), 

by  properly  locating  the  wall  pressure  tap  for  a  given  size 

cylinder. 

Determination  of  thj.s  distance  begins  v/ith  the 

unsteady  Bernoulli  equation 


DV  _  8V   ^9V  ^  _  1  iP  /yN 

Dt  "  at     8s     p  8s  ^'  ^ 


which  reduces  to 


1^  (P  +  pV^/2)  =  p  11^.  (8) 


Upon  integration,  (8)  becomes 


(P  +  pV^/2) 


2 

=  -P  /  11^  ds  (9) 

1     1   ^^ 


or 


2 
?2    +    pV2^/2  -  Pi  -  pV^^/2  =  -p  /  |X  ds.       (10) 


It  is  evident  from  Figure  5  that  Vp  =  0  at  the  stagnation 


point.   Therefore,  (10)  becomes 


2 
^^1-2  "  ^2  ~  ^1  "  "P  •''  It  ^^  "^  pV^/2.         (11) 
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Now,  from  Equation  (6) 


^^1-3  =  ^3  -  ^1  =  '"^   H  <^2) 


and  combining  these  pressure  differences  yields. 


2 

P^  -  P3  =  pvV2  -  p  /  |X  ds  +  pX  II  .         (13) 

2 
In  order  to  satisfy  the  requirement  that  Pp  -  P„  =  pV  /2 , 

it  is  necessary,  considering  (13),  that 


^  5I  =  /'^  ^-  ■  (1^) 


To  evaluate  this  integral,  the  variation  of  V  with  s  along 
the  stagnation  streamline  must  be  known. 

The  complex  potential  for  a  cylinder  in  a  uniform 
stream  is  given  by 


f2  =  V(Z  +  ^)  (15) 


and  the  complex  velocity  is 


2 
W  =  II  =  V(l  -  ^)  .  (16) 


But,  by  definition 


VJ  =  u  -  iv  .  (17) 
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Therefore,  along  the  stagnation  streamline  connecting  the 
points  1  and  2,  one  has 


V  =  0      2 

u  =  V(l  -  ^)  (18) 

s 


since  Z  =  -s.   Substituting  this  result  into  the  integral 
of  Eqn.  (14)  yields 


1  Is 


Carrying  out  the  integration  and  equating  the  results  with 
the  requirement  of  (l4)  yields 


„  8V  _  8V,   .  R^v2  ._-. 


which   simplifies   to. 


2  2 

X  =    (X  +   Y)    +   -q^Vt)    ~   ^  "   V  ^^^^ 


or 


^  +  IxVt)   =   2R  (22) 


where  (X  +  Y)  is  the  distance  to  the  point  2  and  R  is  the 
radius  of  the  cylinder  and  distance  to  the  point  1. 
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Since  the  flow  upstream  from  the  cylinder  extends 
to  relatively  large  distances,  i.e.,  as  X  ->■  °°,  Equation  (22) 
reduces  to 

Y  =  2R.  (23) 

Therefore,  pressure  tap  3  (Fig.  5)   must  be  located  one 
diameter  upstream  from  the  center  of  the  cylinder  in  order 
to  obtain  P^  -  P^  =  PU~/2. 

The  length  of  the  test  section  upstream  from  the 
cylinder  was  sufficiently  large  to  ignore  the  small  influence 
of  the  contraction  of  the  flow  near  the  test  section  on  the 
ambient  flow.   It  is  for  this  reason  that  X  is  taken  to  be 
effectively  infinitely  large. 

5 '      Hydrodynamic  Tunnel  Correction  Factors 

Steady  state  tunnel  correction  factors  as  discussed 
in  Pope  [Ref.  3^]  account  for  the  fact  that  the  bluff  body 
is  not  being  tested  in  an  infinite  medium.   There  is  a 
blockage-effect  correction  to  the  measured  velocity,  which 
in  turn  affects  the  calculated  drag  coefficient.   In  a  usable 
form,  this  correction  factor  is  given  by 


^  =  ri  <f''  '  \  (|>  s  (24) 


where  d  is  the  cylinder  diameter  and  h  is  the  height  of  the 
tunnel  test  section.   This  correction  factor  is  then  used  to 
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modify  the  measured  velocity,  q,  by 


U  =  q(l  +  e) 


(25) 


which  then  modifies  the  uncorrected  drag  coefficient  to 


C^  =  2F/pU  ^d  . 

corr 


(26) 


The  effect  of  this  correction  is  to  increase  the  measured 
velocity  and  thus  to  decrease  the  drag  coefficient. 

However,  with  the  critical  placement  of  the  static 
pressure  tap  3  and  its  proximity  to  the  test  cylinder, 
another  correction  must  be  applied  to  the  measured  velocity 
From  Figure  6,  for  steady  state  flows,  the  ratio  of  local 
fluid  velocity  to  the  free  stream  velocity  is  given  by 


U  = 


2  4 

2  5_.  cos  29  +  ^ 
r  r 


1  - 


(27) 


For  the  geometry  in  this  situation,  the  effect  of  this 
correction  is  to  decrease  the  measured  velocity  and  thus 
to  increase  the  drag  coefficient. 

As  the  effects  of  (25)  and  (2?)  counteract  each  other 
in  the  steady  state  flov;  situation,  and  since  no  information 
is  available  for  time-dependent  flows,  neither  correction 
was  applied.   This  is  particularly  justified  in  the  present 
study  since  the  time-dependent  drag  coefficients  were  compared 
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with  calculated  steady  drag  coefficients,  neither  being 
corrected  for  any  hydrodynamic  constraints. 
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III.   THEORETICAL  CONSIDERATIONS 

Assuming  that  the  unidirectional  periodic  flow  in  which 
the  test  cylinder  is  immersed  is  a  sinusoidal  flow  superposed 
on  a  steady  flow,  the  instantaneous  free  stream  velocity  can 
be  represented  as 


U  =  V  +  V^  sin  2Tr  ft.  (28) 


Therefore,  the  instantaneous  acceleration  becomes 


a  =  ^  =  2iTfV^  cos  2iT  ft  (29) 


and  the  velocity  squared  is 


U^  =  V^  +  2W^  sin  2tt  ft  +  V^^  sin^  27Tft.     (30) 


Now,  in  dealing  with  the  total  force  acting  on  the 
cylinder  in  time-dependent  flows,  it  has  been  most  advan- 
tageous to  consider  it  as  composed  of  a  drag  force  and  an 
inertial  force,  such  as 


F  =  Fj^  +  F^.  (31) 


Then  using  the  concept  of  the  normalized  drag  and  inertia 
resistance  coefficients,  C^.  and  C^,  respectively,  the  total 
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force  can  be  represented  as 


P  =  Cj,!^Apu2  +  Cj^ipA^  §  (32) 


where  A  is  the  effective  cylinder  flow-blockage  area,  A 
is  the  cylinder  cross-sectional  area,  I   is  the  cylinder 
length,  and  p  the  fluid  density.   Upon  substitution  of  (29) 
and  (30)  into  (32),  one  obtains, 


F  =  C^lkApY^   +   ApW^sin27Tft  +  ^ApV^sin^2TTft] 


+  Cj^^[p£A  2TTfV^cos2TTft]. 


2  .  2, 

(33) 


However,  the  results  of  the  works  previously  cited  indicate 
that  a  two  parameter  total-force  equation  does  not  seem  to 
be  capable  of  providing  sufficient  generality  to  cover  the 
conditions  encountered  in  arbitrary  time-dependent  flows. 
Also,  the  past  history  of  the  flow  must  be  considered  along 
with  the  velocity,  rate  of  change  of  velocity,  and  Reynolds 
number.   Therefore,  defining  Cp^(O)  as  the  steady  state 
coefficient  of  drag  v;hich  occurs  at  time  t  =  0  when  U  =  V, 
and  assuming  that  a  normalized  frequency  parameter  X  =  fd/V 
provides  the  necessary  correlation  between  the  resistance 
coefficients  and  the  aforementioned  flow  parameters,  the 
total  force  can  be  written  as 
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P  =  Cp(0)^ApV^ 

+  Cj^(x)[W^Apsin2TTft  +  %ApV^^sin^2TTft  ]    (3^) 
+  Cj^(A)[A^p^27TfV^cos2TTft]. 


The  coefficients  C^CO),  C^(x),  and  Cy,(x)    are  determined 
experimentally.   It  is  obvious  that  the  validity  of  the 
foregoing  analysis  depends  upon  the  existance  of  a  correla- 
tion between  the  said  coefficients  and  the  normalized 
oscillation  parameter,  X,  for  various  values  of  Reynolds 
number,  X,  and  V./V. 


3^ 


IV.   EVALUATION  AND  PRESENTATION  OF  DATA 

A.   RAW  DATA 

As  cited  earlier,  the  differential  pressure  directly 
provides  the  velocity  in  squared  form.   Figure  7  shows  a 
theoretical  curve  of  U^  =  (V  +  V  sin27Tft)^  vjith  V  =  1.5 
and  V.  =  0.5.   The  distances  labeled  p  and  r  occur  at 
t  =  T/^  and  t  =  3T/^,  respectively.   With  these  distances 
knovm,  V,  V.,  and  U  at  time  t  =  T/8  can  be  calculated 
as  follov/s.   Writing  at  the  specific  times. 


t  =  T/4    u^  =  (V  +  V^)^  =  p  (35) 

t  =  3T/4   u^  =  (V  -  V^)^  =  r  (36) 

t  =  T/8    U^  =  (V  +  V^  \J2/2)'^  (37) 


one  has 


V  +  V.  =  >J^  (38) 


and 


V  -  V^  =  \^  (39) 


which  yields 


V  =  (  \/F  +  vy?)/2.  (40) 
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Now,  upon  subtracting  (36)  from  (35),  one  obtains 


P  -  r  =  (V  +  V^)^  -  (V  -  V^)^  =  4W^         (41) 


which  reduces  to 


V^  =  (p  -  r)/4V.  (42) 


Also  upon  expansion  of  equation  (37),  at  time  t  =  T/8,  one 
has 


U^  =  V^  +  \/2W^  +  \1^  (43) 


Figures  8  and  9  display  copies  of  strip  chart  recordings 
of  typical  data  sets  for  low  and  high  frequencies  of  oscilla- 
tion, respectively.   The  upper  traces  represent 

2 
AP2_o  =  pU  /2  v;hile  the  lower  traces  indicate  the  instantan- 
eous total  force..  Due  to  the  fact  that  the  torque  curve  of 
a  butterfly  valve  is  not  symmetrical  about  45°,  nor  the 
discharge  coefficient  (and  thus  the  head  loss)  symmetrical 
about  90°  in  unidirectional  flow  [Ref.  33],  the  velocity 
curve  is  not  a  pure  sinusoid.   Therefore,  the  following 
procedure  was  used  in  reducing  the  raw  data  to  usable  form: 
1.   The  period,  T,  was  measured  between  any  two  corres- 
ponding points  on  the  curve.   Thus  T  was  the  true  period  of 
the  repetitive  signal. 


37 


o 

H 

EH 

< 

tA 
H 
O 
CO 
O 

o 

>-• 
o 

w 

O" 

w 

fin 


o 

EH 

< 

w 
o 

Eh 
O 

o 

H 
P-. 

>H 
EH 


CO 

CD 
U 

:i 
hO 
•H 


38 


o 

M 
Eh 
< 

M 
O 

CO 

o 

fin 
O 

>H 
O 

w 

a 

C5 
M 

Eh 

< 

W 
O 

<c 

Eh 
Eh 

< 
o 

o 

H 

PL4 
>H 
EH 


On 

::! 
bO 
•H 
fin 


39 


2.  The  distances  p  and  r  were  measured  to  the  highest 
and  lowest  points  on  the  pressure  curve,  respectively. 
The  occurrance  of  these  points  were  labeled  T/H   and  3T/4 
irrespective  of  their  distance  on  the  time  axis  since  the 
minimum  and  maximum  velocities  will  occur  at  these  times. 

3.  The  height  at  the  squared-velocity  curve  at  times 
t  =  0  and  t  =  T/8  were  calculated  from  p  and  r  as  if  the 
velocity  were  a  true  sinusoid.   The  corresponding  points  on 
the  curve  were  then  found  at  these  heights,  and  the  base 
line  marked  0  and  T/8. 

^.  Vertical  lines  were  then  drawn  from  the  pressure 
curve  to  the  force  curve  at  the  places  marked  0,  T/8,  and 
T/^,  and  the  height  of  the  force  curve  tabulated  at  these 
points.  Since  the  two  traces  are  made  simultaneously  and 
represent  instantaneous  values,  these  quantities  then  give 
the  forces  which  would  have  occurred  at  times  0,  T/8,  and 
T/^  had  the  velocity  been  a  true  sinusoid. 

5.   The  values  of  velocity  and  force  were  then  calcu- 
lated from  the  calibration  equations.   Reynolds  number 
based  on  the  mean  velocity  v/as  also  calculated.   For  this 
purpose,  an  average  water  temperature  of  70°  F  was  assumed 
for  evaluating  the  density  and  viscosity  of  the  tap  water 
in  the  tunnel. 

B.   EVALUATION  OF  COEFFICIENTS 

As  detailed  in  Part  III,  the  total  force  acting  on  the 
cylinder  is  assumed  to  take  the  form 


^0 


F  =  Cp(0)[%ApV^] 

+  Cj^(X)[W^Apsin27Tft  +  ^ApV^^sin^2Trft  ]     (44) 


+  Cj^(A)[27rA^p£fV^cos27rft] 


After  the  reduction  of  the  raw  data,  at  times  t  =  0,  t  =  T/8 , 
and  t  =  T/4  all  the  variables  with  the  exception  of  Cp^CO), 
C„(X),  and  C^(X)    are  known.   However,  a  system  of  three 
equations  for  the  three  unknowns  can  be  developed.   Re- 
writing (44)  as 

F  =  Cj^(0)K^ 

+    Cp(X)[K2Sin27rft    +   K   sin^2iTft]  (45) 

+   Cj^(X)[K^cos27rft] 

and  evaluating    (45)    at   times    0,    T/4 ,    and  T/8   results   in 

t   =   0  F  =   Cj^(0)K^   +   Cp^(X)K^  (46) 


t    =    T/8      F   =    Cj^(0)K^   +    Cj^(X)[K2    v^/2    +    K^^]    (4?) 
+    Cj^(X)K^    v>^/2 


t    =   T/4      F   =    Cj^(0)K^   +    Cj^(X)[K2    +    K^]  (48) 
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These  three  equations  can  be  solved  simultaneously  for  Cp^(O), 


D 


Cj^(X),  and  C^(\) , 


However,  in  the  range  of  experimentation  of  this 
study,  the  ratio  of  inertia  forces  to  that  of  the  drag 
forces  was  <<  1.   This  fact  led  to  problems  of  instability 
in  the  matrix  used  to  solve  for  the  inertia  coefficient. 
For  this  reason,  Cy.(X)   was  fixed  at  a  chosen  value,  and  the 
values  of  C_.(0)  and  C„(X)  were  solved  using  the  instantaneous 
data  at  times  t  =  0  and  t  =  T/^.   The  two  resulting  equations 
were 


t  =  0       Cj^(O)  =  ^-^^^  (149) 


and 


F  -  C^(0)K, 
T/^    %^^^    =    (K.^-  KJ  (50) 


The  coefficients  C-^(O)  and  C_^(X)  were  calculated  from 
(49)  and  (50)  by  using  two  different  values  for  C^(X) . 
First,  a  theoretical  value  of  C^{X)    =  2  for  unseparated 
potential  flow  and  second,  a  semi-theoretical  and  semi- 
experimental  value  of  Cp,(X)  =  1.5  [Refs.  25,35]  were  used. 
The  resulting  values  of  Cp^(O)  and  C_^(A)  were  essentially 
unchanged.   The  difference  between  C-^(O)  with  C^(X)  =  2 
and  Cpj(O)  v;ith  Cj^(X)  =  1.5  was  negligible  to  the  second 
decimal.   This  also  held  true  for  C^{\) . 
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C.   NORMALIZED  PLOTS 

The  non-dimensional  drag  coefficient  Cp^(O)  is  pre- 
sented in  Figure  10  as  a  function  of  the  Reynolds  number 
based  on  average  velocities. 

The  coefficients  C^(0),  C^(\) ,    and  C„(X)  are  then 
presented  as  functions  of  the  normalized  frequency  parameter, 
X  (Pigs.  11-15). 
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V.   DISCUSSION  OF  RESULTS 

A.   UNIDIRECTIONAL  PERIODIC  MOTION 

Although  the  flow  remained  unidirectional  and  periodic, 
its  oscillatory  motion  was  not  quite  sinusoidal,  as  discussed 
in  Section  IV.   This  aberration  in  the  flovi  appeared  to 
manifest  itself  differently  for  high  and  low  frequencies  of 
oscillation,  irrespective  of  the  value  of  X. 

At  low  frequencies  of  oscillation  (Fig.  8)  the  curve  of 

2 

pU  /2  appears  symmetrical  about  the  vertical  with  the  abso- 
lute values  of  acceleration  and  deceleration  about  equal. 
However,  there  is  a  long  period  in  the  cycle  at  which  the 
valve  appears  to  have  no  effect  on  the  flow.   At  the  highest 
portion  of  the  curve  the  valve  is  in  the  horizontal  position 
with  respect  to  the  conduit  axis  (minimum  head  loss),  while 
at  the  lowest  portion  of  the  curve  the  valve  is  vertical 
(maximum  head  loss).   It  is  thus  apparent  that  the  available 
head  is  not  significantly  changed  when  the  valve  is  about 
±  30  degrees  from  the  horizontal,  causing  a  flat  portion 
in  the  curve. 

At  high  frequencies  of  oscillation  (Fig.  9)j  this  "null 
portion"  of  valve  travel  creates  a  pressure  curve  which  is 
unsymmetrical  about  the  vertical  characterized  by  a  small 
acceleration  and  rapid  deceleration. 

Several  attempts  were  made  to  alleviate  this  problem. 
First,  a  baffle  plate  was  installed  upstream  of  the  valve 
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to  decrease  V.   At  this  time  vertical  ridges  extending 
downstream  were  attached  around  the  holes  in  the  baffle  in 
order  to  decrease  the  magnitude  of  the  deceleration.   Second, 
small  vertical  plates  were  attached  perpendicular  to  the 
valve  faces  so  that  the  valve  might  affect  the  head  loss  in 
the  near  horizontal  position.   Finally,  two  plates  were 
attached  perpendicular  to  the  valve  faces  with  the  same 
diameter  as  the  valve  itself  to  form  a  four-bladed  butterfly 
valve.   Unfortunately,  time  did  not  permit  manufacturing  the 
plates  to  the  same  width  or  contour  as  the  existing  valve. 

It  is  for  this  reason  that  only  those  portions  of  a 
cycle  which  most  nearly  exhibited  a  sinusoidal  velocity 
variation  were  used  in  evaluating  the  various  coefficients. 
A  butterfly  valve  system  consisting  of  several  plates  is 
anticipated  to  alleviate  this  problem  and  yield  a  fairly 
sinusoidal  oscillation  in  velocity  over  the  entire  cycle. 

B.   EXTENSION  AND  COMPARISON  OP  DATA  WITH  THOSE  OBTAINED 
BY  OTHERS 

The  data  obtained  herein  are  presented  together  with 
those  obtained  by  Davenport  [Ref.  36]  at  the  University  of 
Bristol  and  by  Sarpkaya  [Ref.  31]  at  the  Aerodynamische 
Versuchsanstalt ,  Gottingen,  using  a  circular  cylinder  of 
1  inch  diameter  in  an  oscillating  wind  tunnel  within  the 
range  of  Reynolds  numbers  from  5000  to  90,000. 

First,  it  is  apparent  that  the  data  of  Davenport  and 
Sarpkaya  constitute  extensions  of  the  present  data  to  higher 
values  of  the  frequency  parameter,  X. 
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Second,  it  is  clear  that  all  of  the  data  correlate  fairly 
well  with  X  over  a  large  range  of  X  and  Reynolds  numbers 
regardless  of  whether  it  is  obtained  with  air  or  v:ater  in 
a  small  or  large  tunnel. 

Third,  it  is  evident  that  the  effects  of  the  free-stream 
oscillations  on  C„(X)  become  significant  only  when  the 
frequency  of  oscillation  approaches  the  natural  vortex- 
shedding  frequency,  i.e.  when  X  approaches  the  Strouhal 
number  of  about  0.2.   Although  a  theoretical  explanation  of 
the  decrease  of  Cp,(X)  as  X  ->■  St  is  beyond  the  scope  of  the 
present  study,  it  is  not  too  much  of  a  conjecture  to  state 
that  the  formation,  growth,  and  motion  of  the  vortices  are 
significantly  altered  by  the  oscillations  when  the  shedding 
frequency  is  in  the  neighborhood  of  the  oscillatory  frequency 
of  the  free  stream.   When  the  period  of  oscillation  is  several 
orders  of  magnitude  larger  than  the  period  of  vortex  shedding, 
the  flow  behaves  essentially  as  a  juxtaposition  of  instan- 
taneous quasi-steady  states  and  the  total  resistance  can 
be  represented  accurately  by  the  coefficients  C^(0)  and 
Cj^(X). 

Finally,  it  is  evident  that  C  (X)  decreases  sharply  in 
the  range  of  fd/V  from  0.25  and  0.5.   The  reasons  for  this 
are  as  follows.   It  is  a  well-known  fact  that  the  drag 
acting  on  a  cylinder  undergoes  oscillations  with  a  frequency 
twice  the  frequency  of  vortex  shedding.   In  other  words,  as 
the  lift  force  changes  its  direction,  up  or  dov;n,  the  drag 
force  oscillates  with  an  amplitude  of  AC^  about  a  mean  value 
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of  C-^.   Denoting  the  oscillation  frequency  of  drag  with 
f_^,  we  then  have  fi-,d/V  ~  0.44,   The  same  parameter  for  the 
lift  oscillations  is  fjd/V  ~  0.22  v;ithin  the  range  of 
Reynolds  numbers  studied.   Thus  the  data  presented  in  Fig.  l4 
show  that  C^(A)  is,  and  should  be,  affected  in  the  range 
of  fd/V  from  0.25  to  0.5,  i.e.  as  the  drag-oscillation 
frequency  nearly  synchronizes  vjith  the  free-stream  oscillation 
frequency.   Only  a  complete  analytical  investigation,  if 
at  all  possible,  can  explain  as  to  why  C^(X)  should  decrease. 
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VI.   CONCLUSIONS 

The  conclusions  which  can  be  drawn  from  this  study  are: 

1.  The  frequency  parameter  X  correlates  fairly  well  the 
data  obtained  by  various  investigators  regarding  the  resis- 
tance coefficients. 

2.  In  the  range  of  0  <  A  <  0.05,  the  resistance 
coefficients  are  not  materially  affected  by  the  free-stream 
oscillations.   This  range  has  not  been  previously  investi- 
gated.  All  or  most  of  the  previous  studies  have  dealt  with 
oscillations  v;ith  a  zero  mean  velocity,  i.e.  with  X  ->  ». 

3.  The  drag  and  inertia  coefficients  V7hich  are  functions 
of  the  frequency  parameter  X  are  significantly  reduced  in 
the  range  where  the  free  stream  oscillations  nearly  coincide 
with  the  flow  induced  oscillations  of  lift  and  drag. 

h.      The  reasons  leading  to  the  decrease  in  Cy.(X)  and 
C-,(X)  in  the  vicinity  of  X  =  0.4  and  X  =  0.2  respectively, 
can  be  clarified  only  by  a  combination  of  additional  theoret- 
ical and  experimental  investigations  where  the  full  Navier- 
Stokes  equations  are  solved  for  an  oscillatory  flow  about  a 
circular  cylinder  and  v;here  the  characteristics  of  flow  in 
the  near-wake  of  the  cylinder  are  measured. 
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VII.   RECOMMENDATIONS  FOR  FURTHER  STUDIES 

The  following  are  recommendations  for  further  studies 

1.  Improve  the  butterfly  valve  system  so  that  nearly 
sinusoidal  oscillations  may  be  obtained  in  the  range  of 

0  <  X  <  2.0. 

2.  Undertake  a  study  of  the  solution  of  the  Navier- 
Stokes  equations  via  finite-difference  or  finite-element 
methods  for  the  range  of  X  values  cited  above. 

3.  Carry  out  experiments  with  flat  plates  and  other 
shapes  of  bodies  to  ascertain  the  universality  of  the 
parameter  X. 
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3.  ABSTRAC  T 


The  instantaneous  force  acting  on  a  circular  cylinder  immersed  in 
a  unidirectional  periodic  flow  characterized  by  U  =  V  +  V;;^  sin  2TTft 
was  investigated.   The  total  force  was  divided  into  three  components 
consisting  of  a  steady-state  drag;  an  oscillating  drag  which  is  a 
function  of  the  frequency  parameter   X  =  fd/V  ;   and  an  inertial  force 
which  is  also  a  function  of  the  frequency  parameter.   The  three 
resistance  coefficients  associated  with  the  force  components. 


C^(X),  and 


Cj^(X) 


v;ere  calculated. 


Cj^(O), 


It  was  found  that  the  resistance  coefficients  are  not  materially 
affected  by  the  free-stream  oscillations  whenever  the  frequency  of 
the  oscillation  is  lov/er  than  the  natural  frequency  of  the  vortex 
shedding.   Hov/ever,  when  X  approaches  the  Strouhal  frequency  of 
vortex  shedding,  the  inertia  coefficient  Cjvj  is  significantly  reduced, 
and  likewise  when  X  approaches  the  value  of  twice  the  Strouhal  number, 
Cp(X)  is  markedly  reduced.   In  regions  below  this  range,  the  inertial 
force  becom.es  negligible  and  the  total  resistance  could  be 
represented  accurately  by  the  coefficients  CpCO)  and  Cd(X)    with  Cj^ 
fixed  at  its  theoretical  unseparated-flow  value,  2.0. 
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